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Abstract 
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1 Introduction 


It has been known for some time that a D-string ending on a D3-brane looks like 
a magnetic monopole from the point of view of the gange theory living on the D3- 
brane. In ref. [3] it was shown that a fundamental string ending on a D-brane can be 
described as a solitonic solution of the Born-Infeld action of the D-brane. This led to 
some speculation as to whether a D-string ending on a D-brane can be described in a 
similar way. In refs. [4] and [5] it was shown that this is indeed the case for a D-string 
ending on a D3-brane (ref. [5] also contains a discussion of the D5-brane case). This 
D-brane conhguration has been widely studied in the literature. See also, for example 
refs. [6], [7], [8], [9] and [10], and refs. [11] and [12] for related work. 

In ref. [1] we initiated a calculation of the energy radiated during the scattering of 
two D-strings stretched between D3-branes, using the D-strings’ action. Our aim was 
to compare our result with the monopole result, which was shown by Manton and 
Samols in ref. [2] to be 

Erad ~ l.SbmraonVtoc > ~ 1-35 , ( 1 ) 

where Erad is the energy radiated, Etot is the total energy in the system, rrimon is the 
mass of each monopole and V-oo is the asymptotic velocity of each monopole. In ref. [1] 
we discussed the soliton solution to the D-string Born-Infeld action that represents 
two D-strings stretched between two D3-branes. We used this solution to describe D- 
string scattering using the moduli space approximation of ref. [13]. Then we calculated 
equations of motion for perturbations to the moduli space approximation, since the 
perturbations contain the information about the energy radiated. Previous work 
regarding perturbations of the Blon spike can be found in refs. [14], [15], [16] and [17]. 

In this paper we will conclude our investigation with our numerical calculations for 
the energy radiated during D-string scattering. We will seek to solve numerically the 
full equations of motion, rather than using the moduli space approximation. 

The layout of the paper is as follows. In section Owe will review the D3-D1 brane 
conhguration, including the results of our hrst paper, ref. [1]. We will discuss some 
of the implications for our numerical calculations, and make some points about the 
energy of the system. In section Owe will describe the motion of the D-strings in the 
asymptotic limit, when they are far apart. We will discuss how the solutions can be 
split into zero modes and non-zero modes, and we will show that the non-zero modes 
decouple from the zero modes in this limit. It is the energy in the non-zero modes 
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after scattering that represents the energy radiated. In section |4l we will describe our 
numerical calculations to solve the equations of motion, and in section [SI we will use 
these numerical solutions to calculate the energy radiated. Section |S1 contains our 
conclusions. 


2 Describing Magnetic Monopoles Using the Born- 
Infeld Action 

In this section we will briefly review some of the background material, and discuss 
the implications for our calculations 

2.1 Monopoles as Soliton Solutions of the Born-Infeld Action 

We start by reviewing the description of magnetic monopoles as soliton solutions 
of Born-Infeld actions from refs. [4] and [5]. The relevant D-brane configuration is 
D-strings attached to D3-branes. Note that this configuration can be studied either 
using the D3-brane Born-Infeld action, or using the D-strings’ Born-Infeld action. 

Consider the Born-Infeld action for a D3-brane with the magnetic held on the brane, 
Bi, excited, and with a single transverse held, <h, excited. If we look for a static 
solution which minimises the energy, we And that the solution obeys the usual BPS 
equations for a magnetic monopole 


Bi = , ( 2 ) 

where Di is the covariant derivative with respect to the gauge held on the D3-brane 
and i, j = 1,2,3 label the spacelike dimensions of the D3-brane. The simplest solution 
to ((21) is 

= = (3) 

where r is the radial coordinate in the D-brane’s worldvolume. Note that this solution 
for <l>(r) indicates that the D3-brane has been pulled into an infinitely long spike in the 
direction corresponding to the held fl). Calculating the energy and Ramond-Ramond 
charge of this solution, we can conclude that it represents N semi-infinite D-strings 
attached to the D3-brane at the origin. 

On the other hand, this configuration can also be studied using the non-Abelian 
Born-Infeld action for N D-strings. See refs. [18], [19] and [20] for more details of the 
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non-Abelian Born-Infeld action, and ref. [21] for a recent discussion. Exciting three 
transverse scalars in the action, say and and looking for a soliton solution, 

leads to the BPS equations 

+ = (4) 

where a is the D-strings’ spatial direction. The equations (|H) are Nahm’s equations 
from the ADHMN construction of a magnetic monopole (see ref. [22] for a review of 
the ADHMN construction, and refs. [23] and [24] for discussions of how the ADHMN 
boundary conditions apply in this case). The solution corresponding to N semi¬ 
infinite D-strings Tunnelling out’ into a D3-brane is 

*- = ±|. (5) 

where the a* are an iV x representation of the SU{2) algebra. 


We next review the solution to Nahm’s equations which describes two D-strings 
stretched between two D3-branes, which we discussed in our previous paper ref. [1]. 
We defined a new string coordinate, ^ = 2a/L, where L is the distance between the 
D3-branes (note that this coordinate transformation sets the distance between the 
D3-branes to be 2). We introduced the ansatz 
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no summation over i , 


( 6 ) 


where ai are the Pauli matrices. Then Nahm’s equations reduce to 


/(-/2/3 = 0, /'-/3/i = 0, /'-/i/2 = 0, (7) 


where ' denotes differentiation with respect to The appropriate solutions to these 
equations, which were first derived in ref. [26], are 

f If -K(k)dr,{K(m,k) 
sn(K(k)(,k) ■ ■’ sn(K(k)(,k) ’ 

—K(k)cn(K(k)^, k) 
sn(A:(/c)^,/c) 

where K{k) is the complete elliptic integral of the first kind and sn(^, k), cn(.^, k) and 
dn(^, k) are the Jacobian elliptic functions. (See ref. [25] for a review of the properties 
of elliptic functions.) The parameter /c is a modulus with 0 < fc < 1. The fi all have 
poles at .^ = 0 and .^ = 2, which correspond to the D-strings Tunnelling out’ into 
D3-branes. Also, symmetry properties of the Jacobian elliptic functions imply that 
/i and /2 are symmetric about ^ = 1, while /s is antisymmetric. 


/i(e,fc) = 
/3(e,fc) = 
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2.2 Describing Monopole Scattering 

In ref. [1] we explained how to describe the scattering of two D-strings nsing the 
solutions (jHI). 

In the limit k ^ 1 the /* flatten out, and we can write 


/i- K{k) , /2 ~ 0 , /s ~ 0 , as/c^l. 


(9) 


These approximations are accurate, except near = 0 and = 2, where the fi have 
poles. This limit corresponds to the two D-strings being far apart on the xi-axis, 
with positions Xi = ±/i. 

In order to describe the scattering we take the /* with k close to 1 as an initial 
condition. We start the motion by allowing k to depend on time, k —> k{t), such 
that initially the D-strings are moving slowly towards each other. Then k decreases 
with time until it reaches fc = 0, when the configuration is axially symmetric in the 
X 1 -X 2 plane, since /i(^,0) = At this point /i and /2 swap roles, and k 

increases towards 1, so that the D-strings are moving apart along the a; 2 -axis, and 
have therefore scattered at 90°. 

Note that the description of scattering in the previous paragraph relied on the moduli 
space approximation of ref. [13]; the solutions were only allowed to depend on time 
through k{t). In the moduli space approximation the motion of the D-strings follows 
a geodesic in moduli space, so that at any point in time the solution has the form of 
the static soliton solutions (jH)). The D-strings have the same velocity at the end of 
scattering as they did initially, and the potential energy is always constant, so that 
no energy is radiated in this approximation. To calculate the energy radiated we have 
to take into account higher order corrections to the moduli space approximation. We 
can write the full solutions to the equations of motion as 


fdiX) = fid, k(t)) +(id, t) , 


( 10 ) 


where the motion of the D-strings has been split into the zero modes, the /*, and 
the non-zero modes, the e*. The zero modes behave according to the moduli space 
approximation, and the non-zero modes act as small perturbations to the zero modes 
(we can assume that the Cj are small when the D-strings are moving slowly). After 
the D-strings have scattered, the energy in the non-zero modes represents the energy 
radiated during scattering, as we will discuss in section below. 
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2.3 Yang-Mills vs. Born-Infeld 


In order to calculate the energy radiated numerically we will work with the Born- 
Infeld action in the low energy limit, a' —0. There are two ways to take this limit 
(see ref. [1]). The distance between the D3-branes is given by 

L = a'v , (11) 


where, in the D3-brane description, v is the expectation value of the field $, which 
plays the role of the Higgs field. So when we take the limit a' —0 we can either keep 
V fixed, in which case L —0, or we can keep L fixed, in which case v ^ oo. 


The mass of the monopole/D-string is given by 


'^mon TiL , 

9s 


( 12 ) 


where Ti is the tension of the D-string, and gg is the string coupling. Therefore 
the appropriate limit to take is a' —> 0, n fixed, so that the mass of the monopole 
is finite. Unfortunately, taking this limit leads to an action with very complicated 
equations of motion, which we were unable to solve numerically. We will therefore 
take the alternative limit, in which L is fixed, and n —*• oo. The resulting action is a 
Yang-Mills action 

/ oo ^2 

dt / {(fl + + - {^2 - - {^3 - , 

■OO J 0 

(13) 

where Tj = 2ipiai/L. The BPS equations for this action are identical in form to those 
derived from the full action, equation 0, 


- 9 ^ 29^3 = 0 , </92 - ‘^39^1 = 0 , 993 - 991992 = 0 . (14) 

The equations of motion are 

(^1-99"+ 991(992+ 97^) = 0, (15) 

+2 — 9^2 + +2(9^3 + 9^1) = 0 , (16) 

+3 - 9^3 + +3(9^1 + 9^2) = 0 • ( 17 ) 


It is these equations that we have solved numerically, as we will describe in section 21 

Since the mass of the D-strings is infinite in the limit we are taking, the energy 
radiated during scattering will also be infinite, since it is proportional to rrimon- To 
keep our calculations finite we will calculate the ratio of the energy radiated to the 
total energy (which is conserved - see section l2. 4. 2D . since this quantity will not depend 

on TflYnon- 
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2.4 Energy Considerations 

In this section we ontline some general points abont the energy of the Yang-Mills 
system we wish to solve. 

2.4.1 The Energy Densities 

The kinetic and potential energy densities of the system governed by the Yang-Mills 


action (fT!11l are 



K.E. density = 

— (0( + 02 + 03 ) , 

(18) 

P.E. density = 

T 

^ ((v^'l - V 2 V 3 Y + (95( - ^1^3? + {.^3 - V 1 V 2 Y) , 

(19) 


where T is a constant which is determined by the mass of the monopole. We will not 
need an exact expression for T since we will always be dealing with ratios of energies. 

Note that the potential energy density for a solntion obeying the BPS eqnations d 
is zero, as we wonld expect. The potential energy density of a solution therefore 
measures the deviation of the solution away from the BPS solution. 

2.4.2 Energy Conservation 

We will show here that the total energy remains conserved. We consider the Noether 
currents for time translation t ^ t — a, which are 

T 

f = — ((01 + 02 + 0 D + ~ ^ 2 ^ 3 ^ + {^2 ~ ^ 3 ^ iY + {^3 ~ ^ 1 ^ 2 ^) , 

= T((^(0i -f 992952 + 9^303 - • 

The total energy Etot is given by 

Etot = [ fd^ ■ 

Jo 

By current conservation we have 

Etot= , ( 20 ) 

Jo 

which gives 

Etot = T [ 99(01 + 99(02 -f- 99(03 - dt{(pi(p2^3)]l • (21) 


6 


Using equations (l88ll of appendix IB .11 we have 

+ (^9202 + </3303 - = 0(0 

= 0 at ^ = 0 . ( 22 ) 

As discussed below equation (IHKl) in appendix lB.il (pi and p 2 will be evolved in such a 
way that they are symmetric about ^ = 1, and 993 is antisymmetric. The contribution 
to GD at ^ = 2 is therefore also zero. So we have 

Etot = 0 , 

the total energy of the system is conserved. 


3 The D-Strings’ Motion in the Asymptotic Limit 

In this section we discuss the motion of the two D-strings in the limit when they are 
very far apart and moving very slowly. In section EH] we will calculate a series expan¬ 
sion for the position of the D-strings in this limit, which agrees with the monopole 
calculation of ref. [2]. In section EH) we will discuss the decoupling of the D-strings’ 
motion zero modes and non-zero modes in this limit. 


3.1 The Zero-Mode Motion 


When the D-strings are very far apart their interaction is minimal, and we can ne¬ 
glect the contributions from the non-zero modes. The D-strings being very far apart 
corresponds to the limit A; —^ 1 in the solutions ((H]) (see section In this limit we 
can expand these solutions as series in k' = \/l — kP'. We find 




K{k) 

smh.{^K{k)) 


cosh(^A'(fc)) -I- 


1 iK(k) , 

4 sinh(,^A'(A;)) 


^cosh{^ K{k)){Ef 





(23) 


and 

f2{^,k) 


K{k) 

sinh(,^A'(A;)) 


/ l ^K{k)cosh{^K{k)) 

\ 4 sinh(,^A'(A;)) 


(k'f + lcoahHiKimk'? 



( 24 ) 
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and 


Mi,k) 


_jm_( HimcoMcim) . ,,2 

smh.{^K{k)) \ 4 smh.{^K{k)) 


^ cosh‘s{^K{k)){k'Y 


+ 



(25) 


(see Appendix lA.ll for the details of this calculation). The parameter k is unnatural 
to work with because the solutions for /i, /2 and /s depend on it in a highly nonlinear 
fashion. We can use instead K{k), which gives the approximate position of the D- 
strings in the a:i-direction when they are very far apart (see equation inj. This is 
possible in the asymptotic limit because the expansions (I2D - (I2ni) only depend on k 
through K{k). The expansions for /i, /2 and /s in terms of K are best expressed as 
series in We find 




(26) 


and 


f2{^,K) = 


K 


sinh(,^A') 




(27) 


and 


/3(e,^) = 


K 


sinh(,^A') 


1 + 4 


^ a: cosh (^ AT) 
sinh(,^A') 


— cosh^(,^A') ] e 


.-2ir 


(28) 


(see Appendix IA.2I1 . 

In the moduli space approximation the D-strings’ motion is described by allowing 
the modulus K to depend on time. We assume that K is small so that the D- 
strings are moving towards each other slowly (this is necessary for the moduli space 
approximation to be accurate). Energy conservation now takes the form 


T 


dfi 


dt 




df2 


dt 


+ 


df3\ 


(29) 


where we have neglected the contribution of the potential energy because it is very 
small in the asymptotic limit. Here is a constant which represents the initial energy 
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of the D-strings, i.e. their energy when they are an infinite distance apart. Writing 
(EHl) in terms of the parameter K{t), and nsing that /i, /2 and /a only depend on 
time throngh K, we obtain 




dfi 


dK 


df2 


dK 


+ +1^ d( = d^- 


dh 


dK 


Then the expression for K is given by 

k = 

where 


2E 


JiT ’ 


h = 


10 


dK 


+ 


dK 


+ 


dK 


di. 


Using the expansions (jsni) - (jSHI) for /i, /2 and /a in dHS]) we obtain 


I^ = 2[l-- + 0{Ke-^^: 
We differentiate (EDI) to obtain an expression for K 

h 


K = -E 


II ’ 


where 


h = 


dfi\ fd^fi 


dK \dK^ 


+ 


df2\ fd^f2 


dK \dK^ 


+ 


dK 


dK^ 


di 


Again using the expansions JSE!) - (I2HI) in dSSI), we get 

1 


h = 


A'2 


(1 + 0(A'’e-“')) 


(30) 


(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


The leading order terms in the expression 113U give (using equation 113311 1 

K = |u_„| ( 1 - 1 ■ (37) 

where V-^ is the velocity of the D-strings in the asymptotic limit t —> —co. Expanding 
equation ill) as a series in 1 /K, and integrating gives 

K - - hiK + — + O = 'V-oo(t + to) 1 (38) 

where to is a constant parameter corresponding to the freedom to translate the prob¬ 
lem in time. To first order the solution to JSHI) is 


K = V-oo{t + to) . 


(39) 
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We can find higher order solutions by perturbing jHi, substituting the perturbed 
solution back into (El, and solving for the perturbation. The resulting expression 
for K{t) is 


K 


'^—OO 


+0 


+ to) + ^ ln(n_oo 
/ (ln(r;_oot)^\ 

V J ■ 


(t + to)) + 


ln(n-oo(t + to)) 

4v-oo(t + to) 


1 

8v-oo{t + to) 

(40) 


Note that (|40|) is in agreement with the equivalent expression from the three-dimen¬ 
sional monopole calculation, given in equation (8) of ref. [2]. 


3.2 Decoupling of Zero Modes and Non-Zero Modes 

In section we described how the D-strings’ motion can be thought of as being split 
into two parts; the motion of the zero modes, i.e. the motion of the centres of mass 
of the D-strings, and the motion of the non-zero modes, which act as perturbations 
on top of the zero modes. Energy can be transferred between the zero modes and 
the non-zero modes as a result of the interaction between the two D-strings. But 
when they are far apart, and the interaction can be neglected, the zero modes and 
the non-zero modes decouple. We give an argument to show the decoupling between 
zero modes and non-zero modes in appendix IB.1 1 This means that energy can no 
longer be transferred between zero modes and non-zero modes, as we show explicitly 
in appendix IB. 21 It is the energy which has been transferred between zero modes and 
non-zero modes as a result of D-string scattering that represents the energy radiated 
during scattering. 


4 Solving the Equations of Motion Numerically 

In this section we will describe the numerical methods we used to solve the Yang- 
Mills equations of motion, dm) - ini). In section 14.11 we will discuss the numerical 
methods we used. In section lO we will specify the initial conditions we used, and 
in section P~3l we will discuss our boundary conditions. In section we will present 
some graphs of our results for the ipi. 
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4.1 Numerical Methods 


Note from equations JHHD of Appendix IR.1I that near ^ = 0 the take the form 


t) — + o{^). 


( 41 ) 


So the ipi have singnlarities at .^ = 0, bnt the singnlarities are constant in time. In 
order to handle these singularities numerically we removed them by defining the fields 
gi as follows 


9M,t) = + . 


This also implies that 


The equations of motion for the fields gi are 


9i ~ 9i + 9i{92 + 93) 


2 , 2\ (.92 + S'! + ‘^ 9 i 92 + ‘^9193) , (25'i + ‘2g2 + 2 g^) 




92 — 92 + 92 (fi's + 9i) — 


2 , 2\ (fi's + 9i + ‘^9293 + ‘^9291) , (25'i + 2g2 + 25(3) 


93 — 93 + 93(91 + 92) — 


2 , 2^ ( 9 i + 92 ^ 2 g^gi + 25(35(2) (251 + 2^2 + 253) 




+ 

+ 

+ 




= 0 
= 0 
= 0 


(42) 

(43) 

(44) 

(45) 

(46) 


The third and fourth terms in these equations are apparently singular at .^ = 0. 
However, using again the series solntions for the (pi, eqnation (l8^ from appendix 
IR.1 1 we have 


5'i(^)f) — + O(C^) ) 92(^,'t) — bi^ + 0(^^) , 

93(^,'t) = Ci.^ + 0(.^^), 


(47) 


with 

Qji T T Cl = 0 . (4^) 

From this we can deduce that the terms which appear to be singular in (P^ - 
are in fact finite at = 0. 

The (pi also have singnlarities at ^ = 2, but since we will only solve nnmerically for 
the range 0 < .^ < 1, these singnlarities will not affect onr nnmerical calcnlations. 

We solved the equations numerically by evolving an initial configuration in time using 
an RK4 procednre, adapted to two-dimensional partial differential eqnations, snbject 
to certain bonndary conditions (see sections P"^ a, nd P~^ respectively for discnssions of 
the initial conditions and boundary conditions). A detailed discussion of the method 
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we used is given in ref. [28]. Our program was adapted from the RK4 program given 
in ref. [28], and nsed rontines from ref. [29] to calcnlate the Jacobian elliptic fnnctions 
nnmerically. Since the effect we were seeking to observe was so small, it was necessary 
to achieve very accurate results. We did this by using 7 points for the calculation of 
g'l in the RK4 method. We also nsed very small stepsizes; = 0.0001 as the spatial 
stepsize, and dt = 0.00005 as the time stepsize. 


4.2 Initial Conditions 


We start the motion at t = 0 with the monopoles moving tangential to the static 
solutions, so 

= M^,ko) , (49) 

where is chosen so that the D-strings are sufficiently far apart initially. In [2], 
Manton and Samols found that two monopoles cease to interact with one another for 
r > 10. So we take ko = 0.9999999999, for which K{ko) = 12.55264624 /i(^,0). 

At f = 0 the D-strings should be moving towards each other very slowly, so we set 

H^,0) = ko^ , (50) 


where fcg is fixed by the initial velocity of the D-strings as follows. Since fi is ap¬ 
proximately constant in ^ initially, and /2 and /s are approximately zero, /i(l, ko) is 
a good approximation to the initial velocity, Vinu of the D-strings (recall 0 < e < 2, 
so ^ = 1 is the midpoint of the strings). So, having specified Vinu, we can calcnlate 
ko from the following eqnation 


^1(15^0) ko 


dh 

dk 


(l,fco) 


(51) 


The initial conditions for the gi and gi can be dednced from the initial conditions for 
the (fi and (pi respectively using the definitions and iSD. 

Note that the initial confignrations for (pi, (p2 {(fp) and pi, p2 (ps) are symmetric 
(antisymmetric) abont .^ = 1, becanse the fi have these symmetry properties. As 
we discuss in appendix lB.il the (fi will be evolved in time in such a way that these 
symmetries are perserved. We therefore used these symmetries to reduce the numer¬ 
ical computation by solving the eqnations of motion for the range 0 < e < 1, and 
transforming these solutions appropriately about ^ = 1 to obtain solntions for the 
full range of p 
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4.3 Boundary Conditions 

We consider first the left-hand border, ^ = 0. The series expansions (1171) for the Qi 
near ^ = 0 imply 

=^,(0,t) = 0 . (52) 

For points close to the left-hand border it was not possible to use 7 points to calculate 
g'l in the RK4 evolution; we used instead a 3-point calculation for the point next to 
the border, and a 5-point calculation for the next point along. 

To fix the right-hand border we used the symmetry properties of the about ^ = 1. 
These symmetries imply that 

(^;(l,t) = yp'(l,t) = 0, (^ 3 (l,t) = 0. (53) 

These imply for the gi 

= 92i^,t) =-1 , (54) 

g,{l,t) = 1. (55) 

On the right-hand border g^ is therefore fixed, and we have 

g 3 {l,t) = l, ^ 3 (l,t) = 0. (56) 

However, gi{l,t) and g 2 {^,t) vary in time, and so we used an RK4 procedure to 
calculate them. We calculated g'( and g'^ for the RK4 method using the symmetry of 
(pi and ip 2 about ^ = 1, and using 5 points. In the same way we calculated g" and g 2 
for the two points next to the right-hand border. 

4.4 Results for the ipi 

Figures[Tl[2land[Slshow graphs of some of the solutions we obtained from our numerical 
program. 

Figure Q] shows the initial configuration for (pi, p >2 and (^ 3 . From these graphs we can 
see that ipi ~ —77(0.9999999999) = —12.55264624 and 922 ~ 0 and ip^ ~ 0, except 
for the poles at ^ = 0 and ^ = 2 , as we expected from equation ©• 

Figure [21 shows the solutions for ipi, p >2 and ip^, close to the point of scattering. Here 
p)\ ~ p> 2 i which corresponds to the axisymmetric monopole solution (the ‘doughnut’). 

Figure El shows the solutions for p >2 and 923 after scattering, with t = 400. Note 
that after scattering (^1 and (^2 have exchanged roles, as expected. This corresponds 
to the D-strings scattering at 90°. 
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(a) ipi (b) (p2 (c) cp3 


Figure 1: Plots of the numerical solutions for the (ft at t = 0 with Vinu = 0.05. 



(a) ipi (b) (p2 (c) (p3 


Figure 2: Plots of the solutions for the cpi at t = 200 with Vinu = 0.05. 



(a) ipi (b) (p2 (c) 


Figure 3: Plots of the numerical solutions for the ipi at t = 400 with Vinu = 0.05. 


14 



5 Calculating the energy radiated 


In this section we describe the techniqnes we have nsed to calcnlate the energy ra¬ 
diated dnring scattering, nsing the nnmerical solntions for the Qi from the previons 
section. 

The energy densities in terms of the gt are 

K.E. density = ^ (^i + ^2 + qI) , (57) 

RE. density = ^ + g'i + g'i + glgl + glgl + gjgl 

2 

{9i92{gi + 92) + 929?.{92 + 9‘i) + 5'35'i(5'3 + 9i)) 

2 

+ ^(5'i + <72 + S'! + 9i92 + 929z + 9^9i) 

-‘^{9[9293 + 929s9i + 93939i) 

+ -^{91(92 93) + 92(93 + 9i) + 93(91 + 92))^ ■ (58) 

Although there appear to be singularities in the potential energy density dSHD at ^ = 0 , 
all terms are in fact finite when we substitute in the series expansions for the 9i iH) 
(as was the case for the equations of motion (liil) - (lifili b We find 

3T 

P.E. density(^ = 0,t) = ^(®i + ^1 + Cif = 0 . (59) 

In our numerical programs to calculate the energy densities we used four ^-points to 
calculate g[ in dSHl) whenever possible, and we integrated the energy densities using 
Simpson’s rule (see e.g. ref. [30]). 

5.1 Calculating the Energy in the (fi 

In this section we present the results of our energy calculations for the ipi (i.e. the 
total energy in the zero modes and the non-zero modes). We will present all energies 
as ratios to the initial total energy, which we denote 

Eirst we present the ratio of the total energy in the ipi, Etot, to at different times. 
The results are given in table [T] for Vinit = 0.05 and in table El for Vinit = 0.1. If there 
were no numerical inaccuracies in our results this ratio would be 1 at all times because 
the total energy is conserved (see section 12.4.211 . Therefore the order of magnitude 
at which Etot/El'^l^ deviates from 1 at time t gives us an estimate of the numerical 
inaccuracy in our calculation at that time. 
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Time 

77 ' / rninit 

^totj ^tot 

Time 

77 ' / Tpinit 

hjtotj 

0 

1 

400 

0.99999999607 

25 

0.99999999831 

425 

0.99999999663 

50 

0.99999999658 

450 

0.99999999671 

75 

0.99999999639 

475 

0.99999999891 

100 

0.99999999681 

500 

0.99999999925 

125 

0.99999999762 

525 

0.9999999971 

150 

0.99999999706 

550 

0.99999999391 

175 

0.99999999673 

575 

0.99999999648 

200 

0.99999999698 

600 

1.0000000008 

225 

0.99999999709 

625 

0.99999999938 

250 

0.99999999718 

650 

1.0000000012 

275 

0.99999999683 

675 

0.99999999632 

300 

0.9999999968 

700 

0.99999999628 

325 

0.99999999736 

725 

1.0000000017 

350 

0.99999999478 

750 

1.0000000045 

375 

0.99999999561 

775 

1.0000000081 


Table 1: Table showing the ratio of the total energy Etot to the total initial energy 
at different times, with initial velocity Vmit = 0.05. 


Time 

77' / Tpinit 

^totf ^tot 

0 

1 

10 

0.99999999945 

20 

0.99999999947 

30 

0.99999999974 

40 

1.0000000001 

50 

0.99999999987 

60 

0.99999999995 

70 

1.0000000004 

80 

1.0000000003 

90 

1.0000000004 

100 

1.0000000004 

no 

1.0000000004 

120 

1.0000000004 

130 

1.0000000005 

140 

1.0000000006 

150 

1.0000000003 


Time 

TP / TPinit 

^tot/ ^tot 

160 

1.0000000001 

170 

1.0000000007 

180 

1.0000000008 

190 

1.0000000009 

200 

1.0000000008 

210 

1.0000000014 

220 

1.0000000011 

230 

1.0000000007 

240 

1.0000000001 

250 

1.0000000005 

260 

0.99999999947 

270 

0.99999999825 

280 

0.99999999828 

290 

0.99999999844 

300 

0.9999999983 

310 

0.99999999936 


Table 2: Table showing the ratio of the total energy Etot to the total initial energy 
Elol^ at different times, with initial velocity Vinu = 0.1. 
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a) From tabled for Vinu = 0.05, the numerical inaccuracy in the total energy is 
around 10“^. 

b) From tabled for Vinu = 0.1, the numerical inaccuracy in the total energy is 
around 10“®. 

Next we present the potential energy in the ipi in figure |4(aj| for Vinu = 0.05, and in 
figure |4(]^ for Vinit = 0.1. As we pointed out in section 12.4.1 L the potential energy 
measures the deviation of the solution from the solutions to the BPS equations, the 
fi- So the potential energy originates entirely from the non-zero modes. In section 
IB.1 1 we found that the non-zero modes behave like harmonic oscillators when the 
D-strings are far apart. So at late times their kinetic energy is of the same order 
as their potential energy, and so the magnitude of the potential energy in the ipi is 
approximately half the total energy in the non-zero modes. 

a) In the graph in figure |4(a)[ for Vinu = 0.05, we can see that the potential 
energy increases up to around 10“® near the point of scattering t ~ 200. After 
scattering the potential energy decreases back down to order 10“®, which is the 
order of the numerical inaccuracies in this calculation. This suggests that all 
the energy has been transferred back into the zero modes after scattering, and 
therefore no energy has been radiated. 

b) Similarly for Vinit = 0.1, in the graph in figure |4(]^ we find the potential energy 
increases up to the order of lO”"^ around the point of scattering at t ~ 100. Then 
it decreases back down to the order of 10“^ after scattering. Although this is 
slightly higher than the order of numerical inaccuracy, it is still much lower than 
we would expect from the prediction o. which would give Ej-ad/Etot ~ 10 

5.2 Calculating the energy in the non-zero modes directly 

In the previous section we deduced the energy radiated from the potential energy 
of the full numerical solutions for the Qi. In the next section we will calculate the 
energy in the non-zero modes directly. In order to do this we need to separate out 
the non-zero modes Cj from the full solutions (pi, where 

t) = fi{i, k{t)) + ei(^, t) . (60) 
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We will work in the asymptotic limit, when the D-strings are far apart, and we can 
nse K{k{t)) as the modnlar parameter instead of k{t) (see section EH]). We can also 
use the approximations Q. 

In order to separate out the zero modes from the non-zero modes it is necessary to 
calculate K and K for a given numerical solution for the >~pi. A first approximation 
for A', call it K,is K = = 1,^)- This assumes that e 2 {^ = l,t) = 0. To obtain 

a more accurate approximation for K, call it K, we can use the fact that the non-zero 
modes are harmonic oscillators in the asymptotic limit (see section IrTIi . We set 

¥^ 2 ( 15 1) = —K + 5 , (61) 

where 5 = € 2 ( 1 , t) is chosen such that the integral e 2 {^,t)d^ = 0. We then have 

e*(^, t) = t) - fi{^, k) . (62) 

We can calculate an approximation for K and the q using a similar procedure to that 
described above. 

In order to calculate the energy in the e* we have to assume that the zero modes and 
non-zero modes have decoupled from one another, as is the case when the D-strings 
are far apart. Then the kinetic energy density and potential energy density for the 
are given by 

K.E. density = ^(e? + + ^l) , (63) 

RE. density = T+ e'i + e'^) 

+ ^(^1 + ^2 + ^3 + ^1^2 + ^ 2^3 + £3^1) 

+^(^i(4 + 4) + <^ 2(4 + 4) + <^ 3(4 + 4))) ) (64) 

where we have neglected all terms of order and higher in the potential energy 
density 0. 

The graphs in figures El and El show the potential and kinetic energies calculated for 
Vinit = 0.05 and Vinit = 0.1 respectively. 

a) In figure El for Vinit = 0.05, the potential energies are much higher than those 
found in section o for 400 < t < 525 because our calculations were only valid 
in the asymptotic limit. Eor t > 550 the potential energies are of the order 
10“^, which agrees with the results presented in section I01 The we have 
calculated could have originated entirely from numerical errors. 
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Figure 6 



(b) potential energy 

Graphs showing the kinetic and potential energy densities in the e* for 


I 



b) Similarly for Vinu = 0.1, for later times t > 260, the potential energies are of 
the order 10“^, which agrees with the resnlts presented in section I01 


6 Discussion 


We have solved nnmerically the eqnations of motion derived from the Yang-Mills limit 
of the Born-Infeld action for D-strings, starting with two D-strings stretched between 
two D3-branes a long way apart from each other in the D3-brane worldvolume, and 
moving towards each other very slowly. We calcnlated nnmerically the energy radiated 
during the scattering of the two D-strings. 

In discussing our results we should bear in mind that the Born-Infeld action is only 
approximate, and is inaccurate in regions of space which are highly curved. We 
discussed this issue in ref. [1], and gave arguments as to why our results may still be 
relevant. In short, there is evidence to suggest that the solutions to Nahm’s equations 
0 are solutions of the full string theory (see ref. [31] for a recent discussion^). The 
configurations we used were close to solutions of Nahm’s equations, as is demonstrated 
by equation m- It is then reasonable to assume that the motion is accurately 
described by the Born-Infeld action. 

Our numerical results reproduce the 90° scattering which we expected from the com¬ 
parison to monopole scattering. It is encouraging that this follows directly from the 
equations of motion, rather than having to be put in by hand, as we did in ref. [Ij. 

The prediction of ref. [2] that Erad/Etot ~ gives 


■^rad 


E, 


10 


-4 


for V-oo = 0.05 , 


E., 


rad 


tot 


E, 


10 


-3 


for n_oo = 0.1 


( 66 ) 


tot 


(Note that the initial values for our program were Vinu = 0.05 and Vinu = 0.1, rather 
than n_oo = 0.05 and v_oo = 0.1. However, since the D-strings were very far apart 
initially, it is safe to take V-oo = 0.05 and V-oo = 0.1 for our results). Our results for 
the energy radiated are 


^rad 

jpinit 

^tot 


lO" 


for V-oo = 0.05 , 




rad 


Tpinit 

^tot 


10 


-7 


for V-oo = 0.1 


( 66 ) 


(See figures |4(a)] and [Hi for V-oo = 0.05 and |4(b)| and El for V-oo = 0.1). Our results 
therefore suggest that the energy radiated is much smaller then the prediction of 
ref. [2] suggests. Indeed, our values for the energy radiated in (l66l) are approximately 
^We thank A.Gustavsson for bringing this paper to our attention. 
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of the same order as the numerical inaccuracy in our programs (see the discussion 
surrounding tables [T] and [S]). So our results are consistent with there being no energy 
radiated during D-string scattering. It would be nice to be able to support this 
conclusion with further evidence from the Born-Infeld action. 
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A Expansions 


Limit k ^ 1 


In this appendix we derive the series expansions of /i, /2 and /s in the limit k ^ 1 
that we used in section [HI 

A.l Expanding in k'{k) 


In order to obtain /i, /2 and /s as functions of k' we will use the following transfor¬ 


mation of elliptic functions, which can be found in ref. [25]. Under the transformation 


k^k = k'= ^/l-B , 


(67) 


the elliptic functions transform as 


K{k) 

sn(e,k) 


K{k) = K\k) = K{k') , 


dn(e,k) 



cn(^,k) ’ 


1 


( 68 ) 


Substituting into the expressions (jH]) for /i, /2 and /a we obtain 







*^^^^sn(iK(k)e,k') ■ 


(69) 
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We will also use from ref. [25] the expansions for sn(^,k), cn(^,k) and dn(^,k) for 
small k. These are 


sn(e,k) = 


27r /1 , f TT^ \ 5/c^ ■ f \ k"^ . f Stt^ \ 

K{k) Vi \2K{k) J ^ 'el \2K{k) J ^ 64 ) 


cn(^,k) = 


27r /I / TT^ \ 3/c^ f \ k'^ { Svr,^ \ 


dn(^,k) = 


27r 

K{k) 





Substituting these into the expressions for /i, /2 and /s we find 


(70) 


/i(e,fc) 


K{k) 

smh.{^K{k)) 


^cosh.{^K (k)) + 


1 ^K{k) 

4 sinh{^K{k)) 


lcoshKA-(i))(iT 


+ 



(71) 


and 


/2(e,fc) 


K{k) 

smh.{^K{k)) 


/ ui^(fc)cosh(ei^(fc)) 

V 4 sinh(^77(fc)) 


(kT+\cosh^^K{k)){kf 



(72) 


and 

/3(«,fc) 


I<(.k) f, i mk)cosh((K(k)) . y 

sinh(,^i7(/i;)) V 4 sinh(,^i7(fc)) 


jCosh^(^K(k))(k'f 


+ 



(73) 


A.2 Expanding in K(k) 

We can expand K(k) as a series in k', 

K{k) = {\nA-\nk') + ^{hiA-\Tik' -l){k'f+ 0{\Tik\k'Y) , (74) 


which implies 

^-2K(k) ^ itJl + 0{K{k)e-^^^'^^) . (75) 
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Using C3 in m, o and (Hi, we can write /i, / 2 , and /s as series in e to get 

m,K) = 

+ 0 {K‘^e-^^) ) , ( 76 ) 


(cosh{^K) + 4 (- cosh.{^K)\ e 
sinh(^77) V Vsinh(^77) J 


and 


f2{^,K) = 


K 


sinh(,^77) 


V sinh(^77) ' 

+ 0(/C^e-‘“') 


(77) 


and 


f3{^,K) = 


K . /ei^cosh(eir) ^2 


1 + 4 


sinh(^77) \ \ sinh(^77) 

+ 0{K^e-^^) 


— cosh. {^K) e 


-.-2K 


( 78 ) 


B Decoupling in the Asymptotic Limit 

B.l Decoupling of Zero Modes and Non-zero modes 

We consider the D-strings’ motion after scattering, when 

/i ~ 0 , /2 ~-77(t) , /s ~ 0 , (79) 

with K{t) constant. The approximations ( 1791 ) hold trne for all except when ^ is 
very close to 0 or 2, where /i, /2 and /s all contain poles. For now we will work 
with the approximations (Hi; we will consider the effects of the poles later on in this 
section. 

The linearised eqnation of motion for ei from the Yang-Mills action is 

^1 — Cl + ^iif2 + /s) + 2/i(/2e2 + /s^s) = 0 , (80) 

and the equations for 62 and 63 are given by cyclic permutations of (fSn|) . With the 
approximations Hi these equations of motion become 


ei - e'l + K% 

= 0, 

( 81 ) 

^2 — ^2 

= 0, 

( 82 ) 

63 - 6" + Khs 

= 0 . 

( 83 ) 
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Since the D-strings are far apart, K is large, and EJ and (f83l) imply that 


ei = 0 , and 63 = 0 . (84) 

So it seems that in the asymptotic limit the energy in the non-zero modes is entirely 
contained in 62 , which from (m takes the form of a harmonic oscillator, and has 
completely decoupled from the zero mode motion. 

The above analysis is accurate away from the boundaries ^ = 0 and ^ = 2. We now 
consider what happens at the boundaries. Recall that the equations of motion for 
the full fields ^Pi derived from the Yang-Mills action, equation ca , are 

(pi- + + = 0 , Lp2-^2+^2{^l + P>l) = ^ , 

^3 ~ ^3 + + '^ 2 ) = 0 • (85) 

Recall also that our initial condition for the tpi is 0) = /j(^, ko) for some appro¬ 
priate value of ko- The functions fi{^,k), f 2 {^,k) {fsi^yk)) are symmetric (antisym¬ 
metric) about ^ = 1. From the equations of motion (l85ll we can see that the ipi will 
be evolved in such a way that these symmetries are preserved. We therefore need 
only discuss the boundary at = 0 , the same results will follow for the boundary at 
^ = 2 by symmetry. 

Consider the initial conditions for the ipi at ^ = 0. Since (p{^, 0) = fi{^, ko) for some 
ko, the equations of motion give v3j(^, 0) = 0. We also have, at ^ = 0, 0j(O, 0) = 0, 
using the initial condition for the (pi given in equation PI), and that d/i/d/c|(^=o,fc) = 0 
for all k. This means that (pi{0,t) = /^(O, k{t)) at all times because this implies that 
</ 9 j( 0 , t) = (/3j(0, t) = 0 at all times, and so the (pi do not evolve at ^ = 0. The fi have 
the following expansions for small 

/i = -i + oiWe + OK’), /2 = -i + 6iW5 + 0(a . 

/3 — —7 ++ 0(^^) , (86) 

where 

ai(A;) -I- hi{k) + Ci{k) = 0 . (87) 

By continuity there exists a region of small say P, < S, for which the ipi are given 
by 

(pi = —^ + ai{k{t))p,+0{p,^), (p2 =--+ bi{k{t))p, + 0{p,^) , 

(ps = + Ci{k{t))p^ + . 
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( 88 ) 


Since (^i —> as ^ 0, we have that 0 as ^ 0, which gives us the following 

boundary conditions for the e* 


e,(0,t) =e,(2,t) = 0 . (89) 

So we deduce that the boundary conditions are consistent with IHl, and €2 being 
a harmonic oscillator. 

B.2 Energy Decoupling 

We have shown in appendix IR.1 1 that the motion of the zero modes and the motion 
of the non-zero modes decouple when the D-strings are far apart. Therefore we also 
expect the energy in the non-zero modes to decouple from the energy in the zero 
modes in this limit. We will show here explicitly that this is the case. 

The kinetic energy density is given by equation eg. Substituting (pi = fi + we 
find that the coupling between the zero modes and the non-zero modes in the kinetic 
energy is generated by terms like 

2 

Ui di . (90) 

But we have shown in appendix IB. f I that the non-zero Cj behave like harmonic oscil¬ 
lators in the asymptotic limit, and the fi are approximately constant. Then, in this 
limit, 

2 

fiii d^ = 0 , (91) 

and so the kinetic energy does indeed decouple. (The poles of the ft at ^ = 0 and 
f = 2 are fixed, as can be seen from equation JHi, so that fi = 0 at ^ = 0 , 2 . 
Therefore we do not need to worry about the contribution of the poles to jnoD). 

Next consider the potential energy. Substituting ipi = fi + ei into the potential energy 
density, and keeping only terms which are quadratic in e*, we find that the potential 
energy is given by 

T 

RE. = ^ j ((4 - 12^2 - hesf + {cyclic perms.}) d^ . (92) 

Away from the poles the fi are given by the approximations dZHD in the asymptotic 
limit. We have also shown in appendix IB.fl that ei = €3 = 0 in this limit. Then the 
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potential energy is given by 


P-E- ~ ~2 j (^'' 2 )^ + (“/3<^2)^ + (-/ie2)^) Y ^ (4)^'^'^ 

+'^ [ ((4)^ + {- 13 ^ 2 ^ + (-/ie2)^) , (93) 

^ J 2 -S 

where the first and third terms in jnsj take into account the behaviour of the fi near 
the boundaries. Here 5 is a small number, chosen such that the approximations dzi 
are accurate for 5<e<2-5. From the expressions for /i, /2 and /a it can be shown 
that 5 —0 as A; —1. The series expansions dHEI) for the fi imply that 

= -e'(Olc=o + . (94) 


So the contributions to the potential energy from the two boundary terms are given 

by 


3T 


J\2 


+ 


3T 


/ (4)^^e. 

' 2-5 


(95) 


Since these terms are finite, and 5 —> 0 as the D-strings get further apart, the con¬ 
tributions to the potential energy from © are negligible. So the potential energy is 
given by 

r/n p2 

(96) 


T 

RE. = - 


^>\ 2 , 


which has also decoupled from the zero mode motion. 


'0 
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